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ABSTRACT 


Sunshine duration can be estimated using photovoltaic solar cells instead of conventional pyranometers 
or pyrheliometers, which are more expensive and therefore not suitable for low cost measurement 
applications in developing regions. A one-year meteorological dataset from Nicosia (Cyprus) including 
direct irradiance, global irradiance from a pyranometer and global irradiance from a reference PV cell 
was used to calculate sunshine duration following the WMO pyrheliometric method and three 
pyranometric methods by WMO Slob and Monna, Hinssen-Knap and Olivieri. Pyranometric algorithms 
were adapted to the tilted pyranometer and reference solar cell. Main results indicate that all the 
pyranometric algorithms underestimated sunshine duration over the span of a year in Cyprus in 
comparison with the reference pyrheliometric method; and that results between the pyranometer and 
the solar cell were comparable. The PV silicon solar cell is capable of measuring sunshine duration on a 
daily basis with an uncertainty similar to the obtained with a pyranometer when using the Olivieri 
algorithm. 

© 2014 Elsevier Ltd. All rights reserved. 
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1. Introduction 


solar irradiance is above 120 W/m?. This type of measurement is 
often used in low cost renewable energy applications, for example 


The World Meteorological Organisation (WMO) [|1] defines in solar water purification processes such as solar disinfection 
sunshine duration (SD) as the number of hours for which the direct (SODIS) [2-5]. A common approach to SODIS in developing coun- 
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tries uses plastic bottles exposed to the sun to purify water. It 
requires that for a sunny day, 6 h of sunshine is needed to treat the 
water and make it safe to drink. In cloudy conditions, the time 
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required for the water purification increases to 2-3 days. Therefore, 
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a low-cost sensor capable of estimating the sunshine hours would 
be suitable for SODIS water treatment and would improve its 
effectiveness and spread in developing areas [6-9]. 

The aim of this study is to determine if a photovoltaic (PV) silicon 
solar cell can be used to measure sunshine duration, and therefore 
serve as a Sensor for low-cost solar technology applications such as 
solar water purification. To evaluate the suitability of this 
approach, the different algorithms proposed by the WMO for SD 
calculation from direct irradiance, global irradiance and diffuse 
irradiance will be compared using a pyrheliometer, pyranometer 
and silicon cell. 

The algorithms that use only global irradiance from a pyran- 
ometer will be applied to the case of global irradiance data that come 
from a PV silicon solar reference cell [10]. A comparison will be made 
between the sunshine duration calculated by the PV cell with the SD 
values obtained from data from the pyranometer and pyrheliometer. 
The correlation between the data from the cell and the data from the 
pyranometer, as well as possible limitations and possible correction 
factors, will be studied in order to conclude whether the cell is 
suitable or not for the measurement of sunshine duration. 


2. Sunshine duration measurement 


There are different methods to determine the sunshine dura- 
tion according to the WMO [1], including direct measurement 
with the Campbell-Stokes recorder, the pyrheliometric method 
using direct irradiance from a pyrheliometer, or pyranometric 
algorithms using the global irradiance from a pyranometer. There 
are also additional pyranometric methods not adopted yet by the 
WMO but that are well-reviewed in the literature aiming to 
improve the pyranometric algorithm used by the WMO [11-13]. 


2.1. Campbell-Stokes sunshine recorder 


This instrument was introduced in 1880, and is composed of a 
glass sphere that concentrates the solar radiation beam onto a 
graduated paper card that burns according to a sunshine intensity 
threshold. The sunshine duration is read from the total burn length 
[14,15]. The WMO considers that it does not provide accurate data as 
the burns are subjected to errors caused by possible mounting 
adjustments problems and to the fact that the burns depend heavily 
on the card temperature and humidity [1,16]. 


2.2. Pyrheliometric method 


The sunshine duration definition given by the WMO as ‘the 
number of hours for which the direct solar irradiance is above 120 W/ 
m” requires a more accurate method than the Campbell-Stokes 
recorder. In this regard, direct solar irradiance is measured by a 
pyrheliometer held normal to the sun by a sun tracker and 
monitored automatically [17]. A pyrheliometer measures only 
the direct and circumsolar irradiance by using a thermopile with 
a broadband spectral response (typically from below 200 nm to 
4000 nm) and with a narrow aperture. It requires continuous sun 
tracking. The sunshine duration is then obtained by integrating the 
time over the day length during which the direct solar irradiance 
exceeded the threshold of 120W/m?. In summary, the data 
required and the sunshine duration calculation using the pyrhelio- 
metric method are as follows: 


è Data: Direct solar irradiance from a pyrheliometer with a 
resolution of 1 min. 

© Sunshine duration: Period composed by the sub-periods in 
which the direct solar irradiance is above 120 W/m7?. The sub- 
period is 1 min. 


2.3. Pyranometric methods 


Other methods used by the WMO are based on the measure- 
ment of global irradiation using a pyranometer. A pyranometer 
consists of a thermopile with a broadband spectral response; 
similar to a pyrheliometer, but in this case the aperture is not 
narrowed but widened using a semispherical glass dome. If a 
shading ring or a shading ball is used to block the direct radiation 
reaching the pyranometer the diffuse radiation can be obtained. 
The shading ball requires full two-axis sun tracking, whereas the 
shading ring requires weekly elevation adjustment. 

The relationship between direct, global and diffuse solar radia- 
tion is 
I cos 0=G—D (1) 


where I is the direct solar radiation normal to the plane of 
measurement, I cos @ is the horizontal component of the direct 
solar radiation, 0 is the solar zenith angle, G is the global solar 
horizontal radiation, and D is the diffuse solar horizontal radiation. 
If there are two pyranometers available, one for global solar 
radiation and one for diffuse solar radiation then the WMO 
method can be used to calculate the direct solar radiation 
component by using the relationship given in Eq. (1). Sunshine 
duration can then be calculated by filtering for periods where the 
threshold of 120 W/m? is exceeded. Therefore, this method uses 


© Data: Global solar irradiance from a pyranometer and Diffuse 
solar irradiance from a pyranometer with shading ring or 
shading ball and tracker, with 1-min resolution. 

è Sunshine duration: Period composed by the 1-min sub-periods 
in which the direct solar irradiance, calculated as I =(G—D)/ 
cos 0, is above 120 W/m?. 


However, if there is only one pyranometer available, measuring 
global horizontal solar radiation, then the sunshine duration 
calculation is not so straightforward. Several algorithms have been 
proposed by different authors [11-13] that use the global hor- 
izontal and other common parameters, such as the latitude, 
longitude, cloud cover, turbidity, temperature, etc. Some of these 
algorithms are described in the next sub-section. Of these, the 
WMO currently uses the Slob and Monna algorithm. Slob and 
Monna developed this algorithm in 1991 [18]. It uses the mean, 
minimum and maximum of global solar radiation in a 10-minute 
interval. It is based on an estimation of the direct (Eq. (2)) and 
diffuse (Eq. (3)) components for cloudless conditions, which 
depends on the Linke turbidity factor T; [19] (related to the trace 
gases and aerosols in the atmosphere), the solar constant 
(do = 1367 W/m?) and the cosine of the solar zenith angle 
(Mo = cos 0). These estimations are based on a three year dataset 
in the Netherlands (1986-1989) and are as follows: 


I= loexp(— T1 /(0.9 +9.440)) 2) 


where I is the parameterised estimation of direct solar irradiance 
for cloudless conditions, Ig is the solar constant, T; is the turbidity 
factor and jl is the cosine of the solar zenith angle. 


0.2+ 4/3 for 0.1 < pọ < 0.3 
D/G&o= + 03 for fig > 0.3 


where D is the parameterised estimation of diffuse solar irradiance 
for cloudless conditions and Go is the horizontal radiation in the 
atmosphere (Go = loo). 

The algorithm compares the measured global solar irradiance G 
with the lower limit for cloudless conditions, which is Ig +D. This 
comparison is conducted with all the values normalised by Go. 
Fractional values of sunshine f are then calculated for 10-min 
intervals (0 - no sunshine at all, 1 - only sunshine, between 0 and 


(3) 
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Estimation of global and direct radiation for cloudless conditions: 
i (G/G.) gr = H.1/G, E D/G, 
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Fig. 1. The Slob and Monna algorithm to estimate sunshine duration from one-single pyranometer readings of global horizontal radiation. Fractional values of sunshine f are 
calculated for 10-min intervals, comparing with estimated values of direct and diffuse radiations for cloudless conditions [5]. 


Slope = (u;— 1)" bo < 0.3: G/Go< h > f=0 
ee) <GIGg<u, > f=(G/G,- bu, -1)7 
a G/Go= u > f=1 
Bo =0.3: G/Go< h > f=0 
L= G/Go<u > f=(G/Go- hu,- L)! 
G/Go = Uy > f=. 


G/Go 


Fig. 2. The Hinssen-Knap correlation algorithm, showing the linear relationship of sunshine duration with the mean global solar irradiance and the limits established for 
two different intervals depending on the sun elevation angle [2]. 
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1 - partly sunshine, partly clouded), and sunshine duration SD is 
obtained by multiplying f by 10. The complete algorithm is shown 
in Fig. 1. 

In summary, the pyranometric method used by the WMO with 
only a single pyranometer measuring global radiation uses 


è Data: Global horizontal irradiance from a pyranometer. 
© Sunshine duration estimation: Slob and Monna algorithm. 


2.4. Other pyranometric methods 


As we mentioned earlier, there are other pyranometric meth- 
ods developed by researchers at different meteerological agen- 
cies from different countries. One of the most successful and 
accepted algorithms is the Hinssen-Knap algorithm, developed 
by Hinssen and Knap in 2006 by adjusting the Slob algorithm 
[11,20]. The improved algorithm directly relates sunshine dura- 
tion to 10-min mean measurements of global irradiance (Fig. 2). 
There is a lower limit l; for G/Go, and below which there is no 
sunshine (f = 0), and upper limit u;, and above which there is full 
sunshine (f = 1). Between the limits, the sunshine duration is a 
linear function related to the normalised global irradiance. The 
algorithm has two different intervals depending on the sun elevation 
angle ( fy) < 0.3; fg => 0.3). The optimum values for l4, u1, l2 and uz 
need to be established by calculating the pyrheliometric fractional 
values of SD for 10-min intervals and representing against the 
normalised global irradiance G/Go. For the dataset and location 
considered under their study, which corresponds to 1-year data at 
the station of Cabauw (Netherlands, 51.97°N and 4.93°E), the optimum 
values were of l =0.4, u;=0.5, h =0.45 and u= 0.6. These 
values should be calculated for new locations, especially when in 
different climatic areas. 

In a recent report from 2011, Massen |13] has reviewed several 
pyranometric algorithms, including the Olivieri algorithm [21], the 
Slob and Monna, the Hinssen-Knap, the Louche % [22], the 
Campbell [23] and the Glover [24]. He uses the Hinssen-Knap as 
the reference algorithm for comparison, and concludes that 
amongst the other approaches, the most accurate and easy to 
use for calculating the sunshine duration in accordance with the 
WMO definition is the Olivieri one. In 2012, Vuerich et al. [12] also 
presented an evaluation of several pyranometric algorithms. The 
algorithms studied included those by Slob and Monna, and 
Olivieri, among others. They also concluded that the algorithm 
providing the best results, with less uncertainty, was the Olivieri 
algorithm. This algorithm was developed at the Météo France in 
1998 [12,21], and it compares the global irradiance to a threshold 
value that is a function of F, which represents a fraction of the 
global irradiance in clear sky in average conditions of turbidity 
(Fig. 3). The coefficients A and B are specific for each location. 
Météo France calculated an empirical table including the coeffi- 
cients for different location latitudes. 

A description of the main methods to calculate sunshine 
duration using different equipments has been presented: first, by 
means of a Campbell-Stokes recorder; second, using directly a 
pyrheliometer (direct irradiation) and a sun tracker; third, with 
the pyranometric method using two pyranometers (global and 
diffuse irradiations); and finally, using just one pyranometer 
measuring global irradiation, presenting three different algorithms 
in detail. 


3. Estimation of SD from global solar irradiance measured by a 
solar cell 


All the methods to calculate SD described previously are 
suitably accurate but they require expensive equipment such as 


G = pyranometer reading 

h, = horizontal solar angle = 90 -zenithal angle z, 
d = day of the year 

f: sun shines if f= 1, else f = 0 


no sunshine if h, too low 


Gp = 1080*[sin(h,)]" 7° 


; 


F(A,B) = 
A+B*cos(360*D/365) 


G> 
Go*F(a,B) 


Fig. 3. The Olivieri correlation algorithm [4], estimating sunshine duration on 
1-min basis comparison of global horizontal solar radiation with a threshold 
function of a fraction F of the global irradiance in clear sky in average conditions. 
Values of A and B are specific for each location, for this case they correspond to a 
latitude of 44°N. 


a pyrheliometer, a sun tracker or a pyranometer, which is afford- 
able for a meteorological weather station but not for day-to-day 
applications in developing countries. In the case being considered 
in this paper, solar water purification systems in developing 
countries need low cost sensors with an acceptable performance, 
so a trade-off between cost and performance must be achieved. 

A solar cell could be used to measure sunshine duration at low 
cost as it can measure incident solar radiation according to its 
spectral response. However, the solar cell spectral response is not 
broadband but limited according to the energy bandgap (till 
1100 nm for a silicon solar cell). The other limitation is that the 
spectral response is not flat as in a thermopile, but with a 
maximum responsivity in the near-infrared. So the output 
depends on the solar radiation spectrum; and the WMO does 
not advocate the use of sunshine duration detectors based only on 
purely silicon photovoltaic solar cells because these spectral 
variations are a source of error [1]. Other limitations include the 
reduced field-of-view in comparison with a pyranometer and the 
annual degradation of a solar cell (~1% for monocrystalline silicon 
cells) [25-28]. 

In this work, an attempt is made to quantify the effect of these 
spectral variations in comparison with a pyranometer and a 
pyrheliometer when calculating SD. The objective is to determine 
if the SD calculated by a solar cell in relation to a pyranometer is 
well-correlated or not, and if it would be suitable for low-cost 
applications despite the expected performance degradation. 
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3.1. Methodology 


Using a 1-year dataset with direct solar radiation data from a 
pyrheliometer mounted on a sun tracker, global solar radiation 
from a pyranometer and global solar radiation from a calibrated 
silicon photovoltaic solar cell, the sunshine duration was calcu- 
lated using the pyrheliometric method and three of the pyrano- 
metric methods. 

The sunshine duration calculated from the pyrheliometer, 
SDpyrn, Was used as the reference data. The three pyranometric 
algorithms used were the Slob and Monna, the Hinssen-Knap and 
the Olivieri. They were applied to the global data provided by the 
pyranometer, calculating sunshine durations, SDsiop_ pyr, SDHins_Pyr 
and SDoii pyr. They were also applied to the global data from the 
silicon photovoltaic solar cell for comparison, obtaining 
SDstop_sis SDyins_si and SDoii_si- 

As both the pyranometer and the silicon photovoltaic solar cell 
were at a tilted angle facing directly south, and not in a horizontal 
plane, the algorithms were corrected to compare with a tilted 
surface. For Slob and Monna (and so for Hinssen and Knap), the 
estimations of direct normal and diffuse were corrected for a 
tilted surface with an angle / [29] 


Ip =I cos v (4) 

cos y= cos y, cos a; sin P+ sin y, cos p (5) 
1+ cos 

Dg =D ae (6) 

Gp = Ig+ Dg (7) 


where I, is the estimated direct normal at the tilted surface; v is 
the angle of incidence respect to the tilted surface; y, is the sun 
elevation angle; a is the solar azimuth, 2 is the tilt angle; Dz is the 
estimated diffuse radiation at the tilted surface; and Gy is the 
estimated global radiation at the tilted surface. 

The global horizontal extraterrestrial irradiance Gg was also 
substituted for the global extraterrestrial irradiance at tilted surface, 
Gog, using the incident angle v, Go p= lo cos v, for the three 
algorithms. Hinssen correlation optimum limits were established 
for the new solar radiation dataset, obtaining |; =0.1, u = 0.8, 
h =2 and u2 = 0.7. 


3.2. Solar radiation data 


Solar radiation data was obtained from the meteorological 
station and photovoltaic installation from the Photovoltaic Tech- 
nology Group at the University of Cyprus, Nicosia, Cyprus. Latitude 
is 35.2°N and longitude 33.5°E. The direct normal irradiance was 
measured by a Kipp&Zonen CH1 pyrheliometer and the global 
irradiance by a Kipp&Zonen CM21 pyranometer. The calibrated 
photovoltaic solar cell used as a global irradiance sensor is a 
monocrystalline silicon solar cell. Both the pyranometer and the 


Table 1 


calibrated cell are mounted at a tilt angle of 27.5°. Data are 
sampled and stored every minute [30]. 

One-year data were used for this study, from December 2011 to 
November 2012. Data quality was checked in order to first 
eliminate those days with technical problems, such as power 
losses, sun-tracking issues or data acquisition irregularities. Main 
problems leading to wrong or absent recorded data were related to 
the sun-tracker (67 days) and the data acquisition system of the 
photovoltaic solar cell (12 days), usually due to maintenance 
operations and software updates. A second quality control stage 
consisted of filtering to remove solar radiation data that might be 
erroneous, checking the physically possible limits of solar radia- 
tion and the extremely rare limits. 


3.3. Results 


Table 1 presents a summary of the yearly totals of SD for the 
different methods: pyrheliometric, the Slob and Monna pyrano- 
metric algorithm, the Olivieri pyranometric algorithm and the 
Hinssen pyranometric algorithm. The sunshine duration calculated 
by the pyrheliometer is 2171 h. The cumulative difference of the 
pyranometric algorithms over the year is provided, observing that 
the Hinssen and Knap algorithm gives the best estimation, with 
— 145 h(—7%) and — 61 h (— 3%), for the pyranometer and the cell, 
followed by the Slob and Monna algorithm, — 411 h (— 19%) and 
—210 h (— 10%), and the Olivieri algorithm, — 457 h (—21%) and 
—372 h (— 17%). All the algorithms underestimate sunshine dura- 
tion over the span of a year. 

On a daily mean basis, the Hinssen and Knap algorithm gives 
—0.4 + 0.08 h/day (+20%) and —0.17 + 0.11 h/day ( + 65%), the 
Slob and Monna algorithm provides — 1.12 + 0.05 h/day ( + 4%) 
and — 0.57 + 0.08 h/day ( + 14%); and the Olivieri algorithm gives 
—1.25 + 0.06 h/day ( +5%) and —1.02 + 0.07 h/day ( +7%). This 
means that the uncertainty for the Hinssen and Knap on a daily 
basis is extremely high and that either the Slob and Monna 
algorithm or the Olivieri are suitable for measuring sunshine 
duration when using the pyranometer. For the case of the solar 
cell, the most suitable algorithm on a daily basis is the Olivieri. The 
main conclusion of this analysis is that the silicon solar cell is 
capable of measuring sunshine duration on a daily basis with an 
uncertainty similar to the obtained with a pyranometer when 
using the Olivieri algorithm. For the other two algorithms, Hinssen 
and Knap and Slob and Monna, the uncertainty is considerably 
higher than the SD calculation from the pyranometer. 

Fig. 4 shows the calculated daily sunshine duration of the three 
algorithms vs. the sunshine duration calculated by the pyrheli- 
ometer. On the left we can observe the results for the pyranometer 
and on the right for the Si photovoltaic cell. As discussed, the 
Olivieri algorithm gives the better adjustment with the pyrheli- 
ometer for the solar cell. All the algorithms underestimate 
the sunshine duration hours relative to the pyrheliometric sun- 
shine duration. For the case of Slob and Monna, this is already 


Yearly totals of SD for the different methods: pyrheliometric, the Slob and Monna pyranometric algorithm, the Olivieri pyranometric algorithm and the Hinssen pyranometric 
algorithm (h/year), cumulative difference with pyrheliometric SD (h/year) and mean difference (h/day) and standard deviation (h/day). 


Method Instrument SD (h/year) 

Pyrheliometric Pyrheliometer 2171 - 
Slob and Monna Pyranometer 1760 —411 
Hinssen and Knap Pyranometer 2026 —145 
Olivieri Pyranometer 1714 —457 
Slob and Monna PV Si cell 1961 —210 
Hinssen and Knap PV Si cell 2110 —61 
Olivieri PV Si cell 1799 —372 


Difference (h/year) 


Standard error 
of the mean (h/day) 


Mean difference (h/day) 


—1.12 0.05 


—0.4 0.08 
-1.25 0.06 
—0.57 0.08 
—0.17 0.11 
— 1.02 0.07 
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Fig. 4. Daily sunshine duration calculated with the three algorithms (Slob, Hinssen and Olivieri) vs. the sunshine duration calculated by the pyrheliomter for both the 
pyranometer (a, c and e) and the silicon solar cell (b, d and f), showing the correlation between them and the linear fitting. 


comprehensively studied in the literature [11], as this algorithm 
starts considering sunshine when the elevation angle is above 5.7°. 
The Hinssen algorithm lowers this limit to 2.9° and the Olivieri to 
3°. Another reason for underestimation in this particular study is 
due to the tilted surface of the pyranometer, which can result in it 
receiving less light at small elevation angles at sunrise and sunset, 
when the sun can be even behind the pyranometer. This is more 
critical even for the tilted solar cell, as it will not receive solar 
radiation at high azimuth solar angles. However, if the aim is to 
calculate sunshine duration for a particular surface tilted and 


positioned similarly to the silicon solar cell, and with and equivalent 
reduced field-of-view, it will be more accurate to use the solar cell 
than using the pyranometer. This is the same concept as used in 
photovoltaic power plants, using a calibrated solar cell of the same 
technology as the PV modules and in the same position to measure 
solar radiation gives the energy that the photovoltaic modules are 
able to convert into electricity (‘usable energy’), and therefore, 
production estimations and calculations are more accurate [31,32]. 
Fig. 5a gives the frequency distribution for the difference 
between the daily SD calculated by the solar cell using the Olivieri 
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Fig. 5. (a) Absolute frequency of the difference between daily SD calculated with the Olivieri algorithm and the solar cell and the SD calculated with the pyrheliometer (h/ 
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Slob, Hinssen and Olivieri, showing that the Olivieri algorithm gives the better 
adjustment for measuring SD with a photovoltaic solar cell. 


Table 2 

Seasonal pyrheliometric sunshine duration and differences between the sunshine 
durations calculated by the three pyranometric algorithms using the solar cell. 
Autumn months have less total hours due to the reduced number of quality-data 
days due to technical problems. 


Spring Summer Autumm Winter 
SD Pyrheliometric (h) 619 737 272 544 
SD Slob Si cell—SD Pyrh (h) —30 —132 —33 —17 
SD Hinssen Si cell—SD Pyrh (h) 165 25 —119 —133 
SD Olivieri Si cell—SD Pyrh (h) —49 —124 —90 —109 


algorithm and the pyrheliometer. We can see how the SD is 
underestimated as most of the values are below zero. The daily 
mean difference is — 1.02 h and the standard deviation is of 1.4 h. 
Fig. 5b shows the cumulative probability of the daily difference, 
with 95% of the values below +0.25 h of difference. 

Fig. 6 shows the box plot of the differences between the daily 
SD calculated for the three algorithms using the solar cell and the 
SD calculated from the pyrheliometer. It shows again how the 
Olivieri algorithm is the most suitable for the measurement of 
sunshine duration with a photovoltaic solar cell. 

Finally, analysing the sunshine duration calculated by the solar 
cell seasonally, separated into Spring (April-June), Summer (July- 
September), Autumn (October-December) and Winter (January- 
March), it can be observed that the Slob algorithm underperforms 
in the summer months, is similar in spring and autumn and 
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Fig. 7. Seasonal sunshine duration calculated by the pyrheliometer and the three 
pyranometric algorithms using the solar cell: Slob, Hinssen and Olivieri. 


improves in winter (Table 2 and Fig. 7). This agrees with the 
previous studies and analysis [2]. On the other hand, the Hinssen 
algorithm overestimates SD in spring and summer and under- 
estimates SD substantially in autumn and winter, with high 
variation in adjustment. Finally, the Olivieri algorithm under- 
estimates SD over the four seasons, showing greatest underesti- 
mation in the summer and winter months. 

The tilted position of the cell, as discussed earlier, affects the 
performance of different algorithms, as well as the definition of 
the codes for different algorithms. It is important to observe that 
the algorithms were developed mostly in the Netherlands and 
Northern Europe, with different climatic conditions to those from 
the south of Europe, corresponding to the solar radiation data for 
this study. Previous studies have considered datasets with a yearly 
number of sunshine hours of about 1400, and the location in this 
study experiences about 2200 h. It is also expected that turbidity 
values vary considerably from one location to another, so this 
could also affect the performance of the algorithms. 


4. Summary and conclusions 


The objective of this work was to study if a PV silicon solar 
reference cell could be used to measure sunshine duration for low 
cost solar technologies applications. A comparison between the 
standard methods defined by the WMO, using a pyrheliometer to 
measure sunshine duration, and different algorithms when using a 
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single pyranometer, was conducted including the calculation of 
sunshine duration using a solar cell and the pyranometric 
algorithms. 

The evaluation was performed using solar radiation data from 
the meteorological station and photovoltaic installation from the 
Photovoltaic Technology Group at the University of Cyprus, Nico- 
sia, Cyprus, for the period December 2011-November 2012. Three 
pyranometric algorithms were used: the Slob and Monna, the 
Hinssen and Knap and the Olivieri methods. The algorithms were 
adapted to the tilted pyranometer and calibrated photovoltaic 
silicon solar cell from Cyprus. The pyrheliometric method gave an 
annual sunshine duration of 2171 h. The pyranometric methods 
provided annual sunshine durations differences of — 145 h (—7%) 
and —61h (—3%) for the Hinssen and Knap algorithm (pyran- 
ometer and cell); of —411 h (—19%) and — 210 h (—10%) for the 
Slob and Monna; and of —457 h (—21%) and —372 h (—17%) for 
the Olivieri. All the algorithms underestimate sunshine duration 
over the span of a year and the results between the pyranometer 
and the solar cell were comparable. 

On a daily difference mean basis, the Hinssen and Knap 
algorithm had an excessive dispersion and uncertainty in the SD 
values, (—0.4+ 0.08 h/day (+20%) and -—0.17+0.11 h/day 
(+65%)). The Slob and Monna had less uncertainty but still high 
for the solar cell results (—112+0.05 h/day (+4%) for the 
pyranometer and —0.57 + 0.08 h/day ( + 14%) for the cell). Finally, 
the Olivieri algorithm gave a daily mean difference with the 
pyrheliometric method of —1.25+0.06 h/day (+5%) for the 


pyranometer and of — 1.02 + 0.07 h/day (+ 7%) for the solar cell, 
both acceptable results and very similar between them. 

The main conclusion is that the silicon solar cell is capable of 
measuring sunshine duration on a daily basis with an uncertainty 
similar to the obtained with a pyranometer when using the 
Olivieri algorithm. It can measure sunshine duration on a daily 
basis with an uncertainty of 1.4 h/day, which is sufficient for some 
applications. For example, in low-cost solar water purification, the 
device can be overexposed to the sun to reduce this uncertainty. 
Again, this difference and uncertainty value is relative to a 
pyrheliometer, and although it underestimates SD, a silicon cell 
might be more useful as it will give an indication of the real 
sunshine hours that a device with the same characteristics and 
limitations (same position, similar reduced field-of-view) as the 
solar cell is exposed to. 

In terms of cost, the reduction could range from the current 
3000€ per commercial calibrated pyranometer to the cost of a 
silicon solar cell (1-2€/cell plus the cost of calibration, depending 
on the requirements and the lab, up to 100-300€). Another issue is 
that all sunshine duration sensors, whether using pyranometers or 
not, require high accuracy datalogger systems. The current com- 
mercial systems can be up to 3000€, but in another note, low-cost 
dataloggers with high resolution for PV applications have been 
already developed using Arduino [33], with a total cost for the first 
prototypes of only 60€. This cost reduction, along with the good 
performance indicated previously, shows the high potential of 
using silicon solar cells as sunshine duration sensors in developing 
countries. 
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